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FOR: U.S. PROVISIONAL APPLICATION FOR PATENT 



TmJB: ACOUSTIC LOCATION OF A STATIONARY OR MOVING SOURCE 



FIELD OF THE INVENTION 



The pieseixt invention relates to an apparatus and method to locate an acoustic source. More specifically, 
the apparatus and method for acoustic location is based on the method of expanding hemispheres. 

BACKGROUND OF THE INVENTION 

Acoustic location of airborne targets is normally achieved using an analytic solution. For exanq>le, to 
passively monitor the acoustic location of a flying target like the night flight call of a migrant bird, an array 
of sensors (microphones) is used. The signals received firom the sensors are calculated ui time and space. 
The success of this metiiod depends mainly on tiie number of microphones and geometiy of the array of the 
microphones. Often, the accuracy of the results is diminished when a non-standard geometry is used or if 
Hae microphones cannot be adjusted to the same leveL 

One method of calculation is based on time difference of arrival (TDOA). In this case, two sensors will 
receive a signal firom a distant source at difTerent times. This difference in arrival times can be used to 
estimate the source's location (bearing and range). Referring to Figure 1, assuming &r field conditions (Le. 
planar wavefix>nt not spherical) tiie target bearing is determined by Hic delay, Adji^cAtii^x cos9, in fhe 
arrival of the wavefiront at two time-synchronized sensors, Ri and R2. 

The angle G is evaluated by equation (1) 



where At2i is the difference in seconds of the time of arrival of the acoustic signal at the two sensors, Ri and 



6 = cos (cAt2i/x) 



(1) 
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The resulting 9 and flie geometry of the anay define a hyperbolic surface of possible locations of the 
source. The hyperbolic surface is a continuum of points at the calculated angle and distance, in the plane 
perpendicular to the Ime joining R1R2. Anoflier two pairs of receivers are used to determine the hyperbolic 
surfice with respect to each pair. Where ^e three hyperbolic sur&ces intersect, defines the position of 
5 the acoustic source in 3-D, 



The current TDOA calculation method is based on fer field assumptions and hence is inaccurate for a near 
field situation as in the case of acoustic location of a bird in flight In the near field situation, the source is 
assumed to be far enough away fiom the receivers (about 10 times greater than the distance betwera them) 

gl 10 that the acoustic wavefionts can be assumed to be planar rather flian the actual spherical If&rfield 

Q 

y assumptions are valid. At is calculated as shown in Figure 1. Otherwise, a near field situation occurs in 

fU which. At is calculated by acoustic reciprocity starting with spheres located at each microphone and 

m 

^ e3q)andmg these spheres by the speed of sound towards the bird, until they intersect A TniniTrm m of three 

Q 

^ TDOAs are required and therefore four microphones. To survey a sky large enough to find a bird when its 

O 15 location is not predictable, a large beam of acquisition is desked (60^). The radms of the microphone array 
W will determine the height of intersection of the spheres and therefore, the potential sampling heights. A 

Q radius of 50 m provides a sampling height for landing burds 100-300 m. 

H 

It is therefore desired to have an apparatus and method for acoustic detection of a target in flight that can be 
20 acciurate in a near field situation. 



SUMMARY OF THE INVENTION 

It is an object of tiie present mvention to provide a metiiod for determining tiie acoustic location of a target 
in fli^ that is m near a field situation. However, the method is not limitmg to near field situations and 
25 may be used for far field situations, nor is the method limited to movmg target 

Distinguishing features of the method of the present invention mclude tiie incorporation of wavelet analysis 
to allow for the Doppler shift of the entue acoustic signal, and the ability to accommodate non-standard 
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array geometries by using expanding spheres. For examples, non-standard includes uneven microplione 
heights, geometries of more than foxir microphones and witiiout equilateral triangles. 

Conventionally, beamforming is routinely used in antisubmarine war&re to localise a passive acoustic 
5 source. Location is determined by the phase relationships on arrival of a sound wave at an array of 

receivers. Beamforming therefore depends on fiequency. For proper beamforming, the inter-microphone 
spacing, x, determines the maximum frequency of the source,^ that can be resolved according to equation 
(2): 
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where c is the speed of sound. In the location ofbird8,ifone uses c— 340 ms'' and a maximum frequency 
for nii^ flight calls of Inids of 9 kHz, &is method would lunit die inter>miciopIicHie pacing to< 2 cm. 



s 
y 
m 
m 
m 

^ This is undeaiable and can be avoided using tiie time difference ofamval(TDOA) method hi&iscase, 

Q 

B passive localisation of an emitting source involves inteFcepting His source signal at two or more 

s 

Hj 15 microphones. The objective is to determine the source position. 

m 

^ Traditionally, position estimate is achieved using a minimisation procedure using the equations based on 

Figure 2a: 

20 To locate an acoustic source, time difference of arrival (TDOA) or frequency difference of arrival (FDOA) 
can be used. In the case of night flight birds, due to the high frequency of their calls and the speed of sound 
in air, the required distance between receivers does not fit the dimensions of the standard microphones. 
This makes TDOA is tiie more useful method than beamforming. 

25 Using TDOA, the location of the source can be obtained with straight cross-correlation analysis with, for 
example, a night flight call as the matched filter, or by visual inspection of recorded spectrograms. The use 
of wavelet analysis provides a more accurate conelation and thus, more accurate location than either of the 



3 



matched filter or spectrogram analysis approaches. However, the ejcpanding sphere method of die present 
invention can make use of time differences obtained by all three approaches. 



The determination of time differences using wavelet analysis differs from other approaches in that 
correlations are conducted using the matohed filter as well as Doppler-shifted (compressed and stretched) 
fi>rnis of the matched filter. The shifted versions of the matched filter waveform are called baby wavelets 
and a bank or set of these baby wavelets is referred to as a wavelet bank. The recorded tune series is 
compared with the wavelet bank to determine when the signal occurred m each channel. The match with a 
particular baby wavelet is more accurate because the frequency distortion of the received signal can be 
accommodated. Furthermore, by usmg a wavelet bank, the maximum correlation will indicate ^ velocity 
ofthe acoustic target as well as its location. After the match has been fi>und the process is exactly 
con9>aiable to o&er TDOA me&ods ui that the relative time of the peak correlations provides die time 
difference, which is ^plied to the geometry of the microphone array. The present invention provides an 
easier method to deal with complex array geometries. 

More specifically, the method of the present invention is based on the time difference of arrival (TDOA) 
calculation using acoustic reciprocity, i.e. procedure of expanding sphere or hemispheres. The procedure of 
raipanding spheres of the present invention includes the following steps. 

1 . Given a valid (in terms of maximum distance the sound could travel and m^idmiiTn speed of the 
target) set of TDOAs, a mesh of points is constructed on a unit sphere using the following definitions: 




(2) 



(3) 




0 = [l...l...l] ((« + l)xl) matrix 
(4) 

X = cos^*cos^ ((» + l)x (« + 1)) matri 



(5) 
(6) 
(7) 



where 0 and <|> are one dimensional matrices of »+l points, n = 200 (nMnimum, more points for more 
accuracy), ^ increases from the x axis and X, Y and Z are (w+l) by (n+1) matrices comprising the mesh of 
positions on the unit hemisphere (i.e.^ radhis = 1). The elevation angle was restricted to n/2 to n/3 to focus 
the points within the maximum pickup region of the microphones, (note 7i/3 is the lower edge of the 
5 microphone beam for the microphones selected in this experiment This value will depend on the particular 
microphone used and therefore could be more generally stated as the lower edge of the beam of acquisition) 

2. Initialis e tiie size of each sphere based on TDOA values from the correlation analysis: 
XM=D„cX (8) 
10 YM = D„cY (9) 

m ZM = D„cZ (10) 

Q where m is the microphone number (i.e., 1 to 4 in fte test case but tiiis could be increased to any number of 

w 

g7 microphones desired, four is the minimimi required for TDOA), i>m is the delay in arrival of the signal to 

m . 

m each mioxjphone based on the time of arrival at the closest microphone in seconds, c is the speed of sound 

Q IS detennmed by local meteorological conditions in m s'' and X,Y^ were determined by the initial conditions 
^ definmg a mesh on a spherical sur&ce of unit radius (step 1). One is 0, so mitially one of the sets of 

fU equations is identically zero. 

3 

fy 3. Place hemispheres at each microphone position accordiog to the geometry ofthePBFt acoustic 

20 array (i.e., each starting x,y,z position for &e four niicnq>h(nies was positianed on the array): 

(Jfi.yi.2i) = (0,0,0) (11) 

(X2»y2 ,^2) = (0-50,0) (12) 

(pcs^y^ ,^3) = (-50cos(;ir/6),50sm(7r/6),0) (13) 

(X4,y4,z4) = (50cos(^/6),50sin(?F/6),0) (14) 
23 (note the geometry could be changed to any desired arrangement but the test case given here eliminates 

ambiguity over the centre of the array by having a central microphone and three outer microphones situated 

on an equilateral triangle). 

4. Expand hemispheres in time intervals, ti (s) with the growth of each hemisphere determined by c 
30 (m s'') based on local meteorological conditions: 

rsA^ihcit+ti+I^J (15) 
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X„(^,) = /-g^(',)X„ (16) 
Y„(/,) = r^^(f,)Y„ (17) 
^M = r^^{t,)2.„ (18) 
where />„ is the time delay at microphone m in seconds, c is die bcal speed of sound, is die radius at 

die expanded position or next time hiterval at microphone m in metres, x„(tj, yjtif, z„(tj are new positions 
in metres based on die rH.n, expansion at the /, time mterval, where ti = 0.002, and / is the time at the end of 
die previous step in s. 

5. Search die matrix ofpoints on each hemisphere &r intersections: 



m 









M 


M 


M 


•^0+1 


^+1 





(19) 



SI 

m result = {x,y,z\ where XYZ^ r\XYZ^ nXYZ^ r\XYZ^ (20) 

Q 6. Ck)iitinue expanding hemispheres until there is an intersection between all hemispheres within an 



g error region of 1 m\ or the hemispheres have grown beyond a maicinnim ^ » 3 s (note this could be 

Sj 

jy 15 increased if the target was expected to be higher than 1005 m). For the range of values of c at PEPt; Ihis 
^ maximum lumted the search to a radius of about 1000 m (e.g., with c « 335 m s"\ 335 x 3 s = 1005 m). 

m 

The intersection was found using a sub-routine of Matlab® that compares each row (in this case, each 
z) of a pair of matrices and outputs the set of rows common to both. The 'Result" is the intersection between 
20 the wavefionts initiated at microphone I and 2. This is &en compared for overlapping locations with XYZ 
for microphone 3 and 4 and then again for any overlaps between those pairs. 

The main advantage of the above method of expanding sphere is that it enables detection of a target in 
flight located in near field situation. This cannot be achieved with the current method of TDOA 
25 calculation, >^ch is directed to far field situations, and therefore, assumes planar wavefi'onts. 

The acoustic detection system and method of the present invention are ^>plicable for detection of any 
moving or stattonaxy target that emits a sound. This includes, for example, a bkd, a plane, or a snipper. In 
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the latter case, the snipper can be located by the gunshot sound and further tiie bullet's trajectory can be 
determined. 

In one aspect of ^e present invention, tracking and identification of the target in flight can be optimised by 
5 combining acoustic data with radar data. This can be accomplished by fiising acoustic and radar sensors to 
provide a con^>lete system. In the case where the tar^ is a bird, the radar can provide its flight speed, 
biomass and mng-beat characteristics; whereas the acoustic sensor can identify the species of birds. 



Si 10 BRIEF DESCRIPTION OF TH E DRAWINGS 

s 

{jy An embodiment of the present invention will now be described by way of exaniples with refoence to tibe 

m 

f^ following drawings in which: 

m 
.J 

^ Fig^ 1 illustrates the target bearing procedure using time difference of arrival (TDOA) calculation; 

a 

O 15 

m 

fU Figure 2a is a schematic diagram showing the relationship between target and receiver array of 

Q microphones in prior art calculation of a target 

m 

Figure 2 is a schematic diagram of the acoustic system of tiie present invention; 



20 



25 



Figure 3 is a flow chart depicting the method for locating a target based on acoustic data of the present 
imrention; 

Figure 4a is a graph representing a segment of the calibration of time series obtained fiom exaniple 1; 
Figure 4 is a flow chart of the e2q)anding sphere process of the present invention; 
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Figure 5 is a schematic^diagram of a detectioa system of fhe preset that includes a combinatioii of acoustic 
and radar coixq)oneiits; 

Figure 6 is a flow chart depicting the calculation process using radar data to locate a target; 

Figure 7 is a flow chart showing a method for detecting and tracking a target of tiie present invention based 
on acoustic and radar data«; 

Figure 8 is a schematic diagram showing the geometry of &e acoustic array used in experiment described 
in exanqsle 1; 

Figure 9 is a schematic diagram showing the acoustic^radar set-up for the e(;q)6riment described in example 
2. 

DETAILED PESCRIPTION OF THE PREFERRED EMBODIMENT 

The present invention provides an acoustic system and method for acoustic traddng and location of a target 
inflii^ Refeiring to Figure 2, a schematic diagram oftibieacoiustic system is provided. The acoustic 
system includes an array of sensors or microphones connected to one another in a particular geometry. 
Freforabiy, the army has an equilateral shiq[>e wi& a central microphone. For example, using four 
microphones^ a triangular array can be obtained by evenly spacing three microphones at 120^ apart, and one 
microphone centrally positioned This ensures that sound emitted directly over the array centre is captured. 
However, other geometries and additional microphones could be easily incorporated. 

The microphones are furttier connected to a set of pre-anipliflers (not shown), digitiser, and a computer 
where digitised data can be stored. If the target emits a sound while flying above the array of sensors, it is 
captured, digitised and stored in the conqiuter. Once the collection of acoustic data is completed, it is 
analysed to provide the location of the target The analysis uses prior art TDOA wavelet calculation 
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method based on the procedure of expanding spheres of the present uivention (Figure 3). The use of 
wavelet analysis is one novel aspect of the present invention. 



Referring to Figure 3, acoustic field data are analysed as follows. The data is first correlated mth the 
desired range of waveforms such as ballistic or bird wavefonns. Data withm the desired range are farther 
analysed using TDOA The TDOA procedure based on wavelet analysis, locates tiie target signal in each 
microphone channel via peaks in correlation between the received signal and a bank of wavelets which are 
Doppler shifted forms of the target signal at zero Doppler (stationary). Each time a peak correlation occurs 
in multiple channels, the time difference between reception of the same signal provides the data for TDOA 
analysis. These numbers are then submitted to the expanding hemisphere program. 

The expanding sphere process includes the steps of 

1. Given a valid set of TDOAs, a mesh of points is constructed on a unit hemisphere using the 
following definitions: 



where 8 and ^ are one dimensional matrices of n^l points, n » 200 (exceptions givm below), ^ increases 
from the X axis and X, Y and Z are (n+l) by (iH-1) matrices con^rising the mesh of positions on the unit 
hemisphere (i.e., radius » 1). The elevation angle was restricted to to 7c/3 to focus the points within the 
maximum pickup region of the microphones. 




(3) 



(2) 



(4) 
(5) 
(6) 
(7) 



2. 



Initialise Ike size of eac hemisphere based on TDOA values fiom the correlation analysis: 



(8) 
(9) 
(10) 



where m is the microphone nmnber (i.e., 1 to 4), Dm is the delay in arrival of the signal to each microphone 
based on the time of arrival at the closest microphone in seconds, c is the speed of sound determined by 
local meteorological conditions m m and X,Y,Z were detennmed by the initial conditions defining a 
mesh on a spherical sur&:e of nnitxadius (step 1). One Dta is 0, so initially one of the sets of equations is 
5 identically zero. 

3. Place hemispheres at each microphone position according to the geometry of the PEPt acoustic 
array (i.e., each starting x, z position for the four microphones was positioned on the array): 

(0,0,0) (11) 

g 10 (;c„:V2,Z2) = (0-50,0) (12) 

W (x3.:V3,23) = (-50cos(;r/6),50siii(?r/6),0) (13) 

^ (X4,j;4,z^) = (50cos(?r/6),50sin(«r/6),0) (14) 

ru 

Vi 4. i^andhemispheiesm time intervals, /,(s)vvith the grov^ 

S 15 (ms'') based cm local meteorological conditioiis: 
s 

I rs^W=^(^+^.+^») (15) 

nj X„(;,) = r^^(<.)X„ (16) 

^ Y„(<,) = r^^(r,)Y„ (17) 

m Z„(r,) = rj,^(*,)Z. (18) 

20 where is the time delay at microphone m in seconds, c is tiie local speed of sound, n^i ^ ^ radius at 

tiie e:q)anded position or next time uiterval at microphone m in metres, Xt^(Q^ ymM^ ^mM are new positions 

in metres based on die expansion at the ^ tune interval, where /, = 0.002, and ^ is the time at the end of 

the previous step in s. 

25 5. Search the matrix ofpoints on each hemisphere for intersections: 

(19) 

result^{x,y,z\whereX¥Z^r\XyZ^nXYZ^r\XyZ^ (20) 

6. Continue expanding hemispheres until tiiere is an intersection between all hemispheres within an 
30 error region of 1 m^ or the hemispheres have grown beyond a maximum = 3 s. For the range of values of 

10 
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c at PEPt^ this maxinnmi lixoited the search to a radius of about 1000 m (e.g., vnfh c = 335 m s'^ 335 x 3 s 
= 1005 m). 

The intersection was found using a sub-routine of Matlab® that compares each row (in this case, each x,y, 
z) of a pair of matrices and outputs the set of rows common to bofii. The '"result" is the intersection between 
the wave&onts initiated at microphone 1 and 2. This is then compared for overlapping locations with XYZ 
for microphone 3 and 4 and then again for any overlaps between those pairs. 

In one aspect of the present invention, the traddng and location of target in flight can be optimised by 
combining the acoustic system of the present invention with any radar system. Detection by an acoustic 
sensor requires the target to emit a sound signal, and hence combining the two systems enable silent target 
to be detected. Additionally, the range of detection for an acoustic sensor in atmosphere is less than a radar 
sensor. For example, the range of acoustic detection is less than 800 meters for night flight calls of 
nocturnal migrant birds using the microphones selected for the test case; whereas, it is greater than two 
kilometers for radar detection. It is noted tiiat a longer acoustic range is possible but generally at the 
expense of volume of coverage. A larger volume of sky is preferable for targets like birds whose flight 
pattern cannot be predicted. This would also be true for other transient and unexpected targets e.g. 
ballistics. 

In one embodiment of the invention, the acoustic/radar system is used to track birds in flight For the 
purpose of tracking birds, the preferred radar system used is the modified marine radar described in US 
Patent Application 60/304,481, herein incorporated by reference. Briefly, the radar includes a modified 
antexma to enable determination of the bird^s altitude. The antenna is modified by tilting its boresight to 
about 73*' from the horizontal so tiiat the lower half poser pomt of the 26*^ beam is iqyproximately 60'' off 
the ground. The radar is connected to a digitiser and a recorder as shown in Figure 5. Figure 5 also shows 
how the acoustic and radar components are linked. 
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Data collected from the radar are analysed according to the process exemplified in Figure 6. The process 
involves detection, signal processing to form valid tracks and calculation of height of the final 3-D 
positioning along the track. 



5 More particularly, the radar is calibrated according to standard procedures. The radar measurements are 
then converted from polar (r, 9, t) to Cartesian (x, y, z, t) co-ordinates to enable combination thereof with 
those of the acoustic results. The geometric locations are also converted to geographic angles. Tracking is 

j, then performed using standard procedures. The tracks are then validated. For example, if the target is 

moving &ster than) 33 m/s, then the data are discarded ance birds do not typically fry that &st. The valid 

gi 10 track data are then combined with, for example, meteorological data to provide velocity, bearing and height 

J— 7t 

offhetargeL The velocity and bearing calculadomaie obtained using know methods^ 
^ measoiements axe obtain using fhenenlrali^iression method deso^ 

^ 60/304,481. 

a 

fi 

O 15 Briefly, the neutral regression method is based on the assumption that that when the target is directly 

m 

fU overiiead at its closest approach (CPA). 

§ 2W-Rfflm C21) 

20 Where Zest is the estimated target height in meters, R„un is the nrrinfmum measured target range, i.e. CPA of 
the target in meters. Usmg this height estunate, xmdy positions were generated for tiie entire track. When 
this assun^>tion was not valid (i,e., tiie target was not directty above the radar), the position of tiie track at 

was observed to deviate towards tiie radar location (iJoiBet = 0) from the otherwise straight sequence of 
positions. This apparent curving of tracks was noted by Cohen and Williams (1980) when birds were 

25 tracked at the upper edge of a surveillance radar antenna (^ 30** above the radar in the usual horizontal 
orientation). They manually matched the curved track to a computer-generated family of curves. 

If track curvature was observed in this thesis, the elevation angle was decreased in integer iterations from 4» 
= 90^ to 59*" thereby offsetting the radial distance Qience the name, Ro&ed- At each angle, a new track (x, y. 
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z) was conQ>uted until the track was straightened and closer in value to the proper height obtained. When 
the angle was iterated too far, the track was observed to bow out £:om the location of the radar. 

Neutral regression was used to select die track positions that best fit the main axis of tiie trade This fitting 
S technique yields the orientation of the straight line that minimises fiie ooss-linB variance, assuming equal 
uncertainty in jc and (Peaison 1901). 
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The R^sei (range offset (m)) was the value used to correct the altitude and this was calculated £:om ^ as 
follows. 

1 . B^in with the assumption that ibs target is at CPA of the radar - 90^) and iterate fbrougjh values 
of <t> (deg), the elevation angle of the target with respect to the radar» from 90*^ to 59^ The operation 
was applied only at CPA. 

2. Calculate i{o^ values based on values of 

■R^, = (sin(90o) - sm(«>,)) (22) 

3. Calculate the corrected racSus, wilh values of JR^^; 



^offset 




min 



3. Compute x» z at CPA Apply that z to other positions along the track. 

4. Calculate Ri based on z and R^n at each position. 



20 



(23) 
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5. Qdculate the y oo-oxdinates based on J{,: 

jc, cosOg and = sinff^ (24) 

6. Conduct a neutral regression of each set of y points to find the set (defined by that minimises 

the variance from the axis of the track. This value of 4> is annotated (i.e.» ^ optimum). 

7. Ck)nduct a finer sweep using 0.1^ increments in ^ starting at ^qpt*l^ to ^^qpT*'!^, performing a 
regression at each increment until the optimum position is found (minimum length of the semi-mmor axis of 
the principal ellipse). 

The radar and acousdc data are analysed independently to, each, obtain a (x, y, z, t) tracks of tiie target 
The radar data are then cmrelated with the acoustic data to locate overlapping data (Figure 7). The 
overlqyping data confirm the true locations of the target, whereas the non-overlappuig data are discarded. 
This does not mean the independent locations are incorrect It means only that die target did not provide an 
acoustic signal being tracked by the radar or a radar return while being detected acoustically. 

Data within the desired range of waveforms are then analysed using TDOA. If the TDOAs do not make 
sense given die irig xi n^^T n^ distance between microphones and the maxiimnn speed of the target then 
discard, o&erwise the data is fiirther analysed using Ifae method of expanding sphere of the present 
invention (equations 2 to 20) to obtain an (x, y, z, t) location of the target 

The s^or fiision is not limited to radar. Rather, the acoustic method can provide target identification for 
other sensors such as infirared. 

EXAMPLE! 

Acoustic Detection of a Known Moving Source 
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A cyclist lode past an acoustic array while continuously blowing a whistle. The array is shown in Figure 8. 
It has the source in one quadrant of the array, against a surveyed 20. 13 m height above the array and R is 
the receiver. 

Acoustic data collected from the sounds emitted by the whistle were analysed according to the method set 
out in Figure 3. This includes, namely, correlating the field data with whistle wavefomos, thus filtering 
undesired sounds. The zero-Doppler waveform was recorded from a stationary source blowing into central 
microphone. This signal was then digitised and isolated. The whistle wavelet bank was then created by a 
sine frinction compression and expansion of the zero-Doppler (stationary) waveform of the whistle signal. 

The data were then analysed using the time difference of anival (TDOA) method based on wavelet 
analysis. Data that were not possible given the geometry of the microphone array and the speed of the 
target were discarded, and tiiose that were within the expected speed of sound under the conditions of the 
experiment {tenq)erature, wind, target flight behaviour) were further analysed using the method of 
expanding spheres of the present invention (equations 2-20, Figure 4). 

By graphmg the magnitude of the correlation between the signal and wavelet bank, on a plot of the scale of 
tiie Doppler-shifred wavelets versus time, the slope of high correlation ridges indicated the direction the 
target was moving relative to tiie noicrophones as shown in Figure 4a. Figure 4a represents a segment of 
the calibration time series when a cyclist passing by the acoustic array blew a whistie (f(p3kHz, sanapling 
rate=44kHz). The whistle arrived at Channel 1, 3, 2 and 4. The arrows mark the positions used to calculate 
TDOA. Magnitude of the convolution is represented by colour brightness with yellow>red>blue. 

When the target approadied the microphone, tiie slope was negative. The slope was infinite as tiie target 
25 went overtop of the miax)phone, and negative as the target passed tiie microphone. This provides a vaore 

accurate TDOA for location of the target than strai^t matched filter correlation analysis. The location with 
the matched filter would be difficult as the proper correlation peak could be any yellow along the line 
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through scale » 1.0. Note especially, that the coiielation value at CPA (indicated by the axiow) would be 
misleadingly low (less intense in colour). 

The scale of the wavelet which was maximally correlated with the signal, gave a measure of the velocity of 
the target relative to the receiver. In die test case using a cyclist travelling at a known speed, the wavelet 
analysis confirmed this speed (Le., scale conversions gave relative velocities of -5.8 m s*\ -7.8 m s'^ and 
13.8 m s'^ The cyclist's reported speedometer reading was 5.6 to 13.9 m s''). 

The scale provides further inforaiation on the target position (up or down Doppler) relative to each 
microphone (negative velocities so approaching microphones 1, 2 and 4; and positive velocity so leaving 
microphone 3), and therefore, direction of movement The direction of the cyclist was also confirmed by 
segmenting the recorded time series and calculating the change in position over one second interval of &e 
ten seconds time series. The cyclist was correctly determined to be moving east-south-east (110^) (Figure 
8). 

EXAMPLE! 

Detection of Birds Using Acoustic/Radar System 

Swainson's thrush (SWTH) and zeep noctural migrant birds were used to demonstrate the capability of the 
acoustic/radar system. 

The experiment was conducted at Prince Edward Point (PEPt), Picton Country (43**56'00"N, 76^52*00*'W, 
43''45'30"N, l&'S 1 *45*'W) because noctural migrants are known to land there in large xmmbers. The site is 
relatively open, quiet, protected from human interference, and provided landing habitat for noctural 
migrants. The acoustic/radar set-up is represented in Figure 9. 

San^ling was divided into three tune periods: dusk, midnight and dawn. Target behaviour at dawn was 
expected to be significantly different from the other two time periods. At dawn, the flight of landing birds 
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were expected to be reduced in altitude and opposite the mean flight direction of birds sampled at midnight. 
Birds at dusk were expected to climb in height and fly grouped in the same general direction as midnight 
taigets. As the exact timing of the onset of landing was not known, sampling began one hour before dawn. 
Radar and acoustic samples were collected simultaneously for one hour for all diree times periods. 

The collected radar data were analysed according to the metibod shown in Figure 6 to provide tiie results in 
Table 1. 
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Table 1: Locations of SWTH and ZEBP Targets Based on Radar Data 



Date 


Period 


Returns^ 


Intensity^ 


lnt„SD 


Speed* 


Sp^SD 


Altitude'* 


Aft_SD 


n 


Direction® 






28-Aug 


Dawn 


4.3 


195 


17J2 


10 


3.6 


237 


57.9 


-T9' 


SE 


2.254 


0.979 




Dusk 


3.6 


222 


17.4 


21 


6.1 


256 


143.8 


19 


SE 






29-Aug 


Midnight 


4.1 


216 


19.5 


14 


2.9 


227 


88.3 


17 


SB 






Dawn 


4.2 


199 


15.2 


20 


1Z7 


270 


78.3 


18 


SW 


Z413 


0.941 




Dusk 


3.5 


210 


19.2 


18 


8.3 


258 


122.9 


34 


SE 


u% 




12-Sep 


Midnight 


3.1 


238 


11.7 


19 


9.5 


188 


36.3 


41 


SW 






Dawn 


6.3 


213 


18.1 


16 


8.2 


229 


61.3 


39 


SW 


2.28 


0.917 




3.6 


222 


18.6 


19 


8 


277 


79.7 


49 


SW 


2.548 


0.774 






, 3.8 


172 


14.3 


18 


8.4 ■ 


193 


52.4 


41 


NW 


2.69T 


0.884 






3.9 


210 


16.8 


17 


7.5 


"237 


Tar 


31 


■"sW 


3.768 


0.993 



^ average number of pulse returns per track 

^ average intensity of radar pulse return (±SD in next column) 

average (±SD in next column) of all tracks during the sampling period and date in m s*' 
^ average (±SD in next colunm) of all tracks during the sampling period and date in m 
^ quadrant including most track directions (NE,SE,SW or NW) 
^ Shannon diversity index (fis) (Appendix C.7» Table C-4) 
B Shannon's equitability (Appendix a?. Table C-4) 



The acoustic data were analysed according to the mediod shown in Figure 3 to provide the results in Table 
2. 
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Table 2: Locations of SWTH and ZEEP Targets Based on Acoustic Data 



Date 


Time of nfc ^ 


Species 


x(m) 


y{m) 


z(m) 


\lcA fm/^^ ^ 


Range (m) 


/\zimutn {oeg) 


29 Aug 


5:02:18:227 


zeep 


-104 


-99.6 


100 


10.6 


175 


136 




5:05:49:806 


zeep 


■^48 


-14 


438 


33,7 


559 


178 




5:05:49:993 


zeep 


-25.6 


-3.6 


690 


33.7 


692 


172 




5:11:37:085 


swth 


-71.6 


-41 .6 


lO'ie 






10U 




5:11:51:309 


zeep 


-65 


37.5 


164 


22.7 


97 


2o0 




5:12:59:503 


swth 


-36,6 


0.4 


366 


Q 7 


ODO 


lO 1 




5:14:18:638 


swth 


48.4 


-29.6 


750 


lO.f 


r 






5:14:20:196 


swtn 




-O.O 


/ O 




107 


177 




5:14:43:733 


zeep 


-9.6 


6.4 


146 


Q Q 
0.9 




Z I** 




5:14:43:943 


zeep 


-35.6 


16.4 


4 il Q 




lOO 






5:35:37:01 1 


swth 


-/ 


o 




WW 


7fi 


1*9 


18 Sep" 


5:02:08:409 


swth 


-18 


12 


79 


36.4 


91 


191 




5:02:10:775 


swth 


-46 


28 


81 


33.7 


111 


199 




6:02:11:227 


swth 


-44 


24 


96 


36.4 


120 


196 




5:08:34:349 


swth 


44 


-26 


70 


14.2 


78 


58 




5:08:35:248 


zeep 


-43 


-26 


65 


14.2 


99 


157 




d:HS:Sfl:W2 


swth 


26 


-10 


100 


16.5 


101 


88 




5:19:04:308 


zeep 


-268 


236 


220 


25.9 


439 


218 




5:19:05:947 


swth 


6 


47 


124 


10.8 


192 


246 




5:32:33:460 


swth 


101 


-6 


103 


12.0 


^28 


7 




5:39:27:226 


swth"- 


9 


-4 


94 


16.5 


96 


155 




5:39:51:112 


swth 


20 


^2 


81 


^8.> 


98 


96 




5:39:51:644 


swth 


16 


-24 


76 


14.2 


81 


A 

109 




5:39:52:825 


zeep 


10 


-6 


68 


10.6 


70 


148 




5:39:55:197 


swth 


50 


49 


78 


21.1 


94 


298 




5:41:31:123 


swth 


-26 


7 


105 


23.5 


117 


184 




5:41:34.784 


zeep 


19 


-23 


155 


28.5 


157 


104 




5:41:37:557 


zeep 


-28 


61 


71 


28,5 


106 


227 




6:41:37:901 


swth 


-23 


-14 


69 


16.5 


88 


160 




5:42:24:331 


zeep 


38 


-1 


81 


25.9 


82 


20 




5:42:25:387 


zeep 


0 


-27 


185 


23.5 


188 


130 




5:42:26:113 


zeep 


0 


47 








240 




5:42:26:582 


swth 


-15 


-9 


60 


9.9 


74 


163 




5:42:26:605 


swth 


0 


192 


141 


18.4 


237 


262 




6:42:33:561 


swth 


0 


-1 


66 


0.0 


70 


170 




5:42:35:274 


swth 


-35 


-21 


67 


10.6 


94 


158 




5:42:35:367 


swth 


0 


0 


69 


3.7 


74 


172 




5:42:35:749 


zeep ^ 


0 


-1 


68 


25.9 


73 


170 




5:42:37:208 


swth 


25 


14 


79 


5.4 


8G 


267 




5:57:09:902 


zeep 


20 


6 


78 


14.2 


78 


203 



same species 



■ both species 

^Tlme of night flight call (hour.mln:sec:msec) ^ proportion zeep = 55% 
^ Maximum absolute value " proportion SWTH = 62% 
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The two sets of data were tben correlated to locate overl^ping data. The overlapping data confirms 
locations of the targets (Table 3). 



Table 3 : Locations of SWTH and ZBBP Targets Based on Acoustic/Radar Data Fusion 



Date 


Period 


Parameter 


Radar Track Positions 


Acoustic Ijoc 




Data Fusion 


1 


2 


3 


4 


1 5 


1 




Speoes 






Ht (m) 


12-Sep 






505 30 


5:05-33 


5.05*35 


5:0537 




5-05-50 


50 - 38 = 14 


zeep 


11 


254 


435 


Dawn 


Rsnge 


413 


435 


435 


443 


435 








73 


83 


88 


95 








12-S^ 




Tim© 


534:29 


5"34-32 


634.34 


'-•^^^ 




5:34.37 


37 - 27 = 10 


ewth 


10 


45 


488 


Dawn 


Range 


510 


495 


468 






640 






Aziniutti 


165 


162 


160 






1 /z 




18-Sep 


Dawn 


Time 


5.13:69 


5-14.01 


6:14.04 


5:14.08 




6.14:08 


d-0s8 


swth 


15 


225 


98 




180 * 


165 


150 


160 


m. 


140 




MZImUUI 


358 


340 


320 


303 




OOo 




18-8ep 


Daiivn 


Time 


517:54 


5.17-^ 


5:17:58 


5:18*01 


5:18 03 


5:18 04 


64-5S°9 


zeep 


12 


252 


218 


Range' 


218 * 


225 


240 


255 


285 


203 




AZiinuin 


210 


215 


220 


223 


225 


214 




18-Sep 


Dawn 


Time 


5:40:55 


5:40-56 


5:41*02 






5*40:26 


28 -59 = -33 


swth 


7 


292 


83 


Range^ 


150 


157 " 


* 195 






188 




Asimuth 










262 




19-Sep 


Dawn 


Time 


5:07:06 


6.07:10 


5:07:13 


5:07:15 




507 09 


9-15a-7 


swth 


37 


120 


161 


Ranpe* 


186 


285 


315 


435 * 




473 




Azimuth" 


55 


63 


85 


93 




92 




19^p 


Dawn 


Time 


5:11-02 


5-11:05 


5.11:07 






5:11:13 


13-3=10 


thrush 


18 


352 


165 


Range' 


165 * 


^ 158 


173 


M 




188 




Azimuth" 


60 


70 


113 






63 




l&^p 


Dawn 


Time 


511:07 


5:11.09 


6:11:12 


5:11:14 




5:11:13 


13-10 = 3 


thrush 


25 


167 


135 


Range" 


* 188 


180 


173 


165 




270 




Azimuth" 


225 


233 


238 


285 


!^ 


214 






19-Sep 


Dawn 


Time 


5.11:16 


5-11:19 


511:21 


511-.23 




5:11.23 


23-18«5 




13 


347 


143 


Range"* 


158 


* 158 


150 


143 




90 




swth 


Azimuth" 


100 


76 


68 


65 




84 






19-Sep 


Dawn 


Time 


511:26 


5:11:28 


5:11S0 






5-11-24 


24-23 »3 




26 


62 


293 


Range' 


* 293 


300 


308 






360 




swth 


Aarouth" 


330 


340 


35D 






313 






19-Sep 


Dawn 


Time 


5:28:43 


5:26-50 


»»=^«.«^ 






5-26.43 


43.47s^ 




15 


318 


113 


Range^ 


153 


* 153 








143 




thnish 


Azimuth" 


75 


20 








36 





'inm<1bln»7.5m} 
*'Inde8(1binBO5de£0 

^ time of acousbc location minus time on radar track In seconds 
* esOmated timing of acoustic signal 



Tables 1, 2 and 3 show a variety of birds, which demonstrate that fbs invention works for a large range of 
birdies. 

While the foregoing embodiment of the present invention have been described and shown, it is understood 
that all alternatives and modifications may be made and fidl within the scope of the invention. 
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ABSTRACT 

A method of acoustic location of a stationary or moving source is disclosed. The method is based on the 
time difference of arrival (TDOA) calculation using acoustic reciprocity, i.e. procedure of expanding 
spheres or hemispheres. A detection system is also disclosed, which includes the combination of acoustic 
and radar detection ^paratus. 
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